Early anisotropic hydrodynamics and the RHIC early-thermalization and HBT puzzles 
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We address the problem if the early thermalization and HBT puzzles in relativistic heavy-ion 
collisions may be solved by the assumption that the early dynamics of the produced matter is 
locally anisotropic. The hybrid model describing the purely transverse hydrodynamic evolution 
followed by the perfect-fluid hydrodynamic stage is constructed. The transition from the transverse 
to perfect-fluid hydrodynamics is described by the Landau matching conditions applied at a fixed 
proper time r tr . The global fit to the RHIC data reproduces the soft hadronic observables (the pion, 
kaon, and the proton spectra, the pion and kaon elliptic flow, and the pion HBT radii) with the 
accuracy of about 20%. These results indicate that the assumption of the very fast thermalization 
may be relaxed. In addition, the presented model suggests that a large part of the inconsistencies 
between the theoretical and experimental HBT results may be removed. 
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I. INTRODUCTION 

The experimental results obtained at the Relativis- 
tic Heavy-Ion Collider (RHIC) are nowadays interpreted 
as the evidence that the matter produced in relativis- 
tic heavy-ion collisions equilibrates very fast (presumably 
within a small fraction of 1 fm/c) and its behavior is very 
well described by the perfect-fluid hydrodynamics [l|-[3]- 
The fast equilibration and perfect-fluidity are naturally 
explained by the assumption that the produced matter is 
a strongly coupled quark-gluon plasma (sQGP) @. An- 
other possible explanation assumes that the plasma is 
weakly interacting, however the plasma instabilities lead 
to the fast isotropization of matter, which in turn helps 
to achieve equilibration [10(. 

The concept of a very fast equilibration is difficult to 
reconcile with the models of the early stages of the rela- 
tivistic heavy-ion collisions, which typically use the ideas 
of color strings or color-flux tubes. The system produced 
by strings is highly anisotropic; the pressure in the direc- 
tion transverse to the collision axis is usually much larger 
than the longitudinal pressure. 

Interestingly, the decays of fluctuating strings produce 
matter that looks as if it was thermalized only in the 
transverse direction, i.e., the transverse-momentum spec- 
tra of partons are exponential [ll[ . This feature inspired 
the construction of the hydrodynamic model where only 
transverse degrees of freedom are thermalized, while the 
longitudinal motion is described by the free streaming 
[l2| . We shall call this framework the transverse hydrody- 
namics. The first calculations using this idea showed that 
it was consistent with the data describing the transverse- 
momentum spectra and the elliptic flow coefficient V2 of 
the pions measured at RHIC [ID, [l3|. Further studies, 
where the transverse hydrodynamics was combined with 
the statistical hadronization model, showed that one may 
reproduce in this framework also the HBT pion radii in 



a quite successful way [3, EH . 

The concept of transverse hydrodynamics allows us to 
avoid the problem of early thermalization. Neverthe- 
less, the RHIC data suggests that the matter is even- 
tually fully thermalized (by this we mean that the local 
thermodynamic equilibrium with the locally isotropic mo- 
mentum distribution is achieved). Thus, we expect that 
during the evolution of the system the early transverse- 
hydrodynamics stage should change into the stage de- 
scribed by the perfect-fluid hydrodynamics (with possi- 
ble small corrections from viscosity and other dissipativc 
effects). In this paper we present the hybrid model that 
incorporates such a change. In this way we generalize the 
methods and concepts introduced earlier in Refs. [L214I1 ] 

We assume that the transverse-hydrodynamics stage is 
valid in the proper time interval t\ < r < Tt r . The proper 
time Ti defines the initialization time for transverse hy- 
drodynamics, and T tr is the transition time from the 
transverse hydrodynamics to the perfect-fluid hydrody- 
namics 1 . This transition is described with the help of the 
Landau matching conditions. Both t; and r tr are the pa- 
rameters of the model. For r > r tr the system is treated 
as the perfect fluid. The freeze-out takes place on the hy- 
persurface of constant temperature, and the Monte-Carlo 
method is used to calculate the soft hadronic observables: 
the transverse-momentum spectra of pions, kaons, and 
protons, the elliptic flow of pions+kaons and protons, 
and the HBT radii of pions. 

As stated above, the idea of transverse hydrodynamics 
helps to solve the problem of very fast thermalization. 
There exists, however, another important problem in the 
interpretation of the RHIC data, namely, the HBT puz- 



1 We shall refer to this transition below shortly as to the 
2D — > 3D transition. 
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zle. This name refers to general problems connected with 
the correct reproduction of the HBT radii in the hydro- 
dynamic models. One possible solution of this problem 
suggests the use of the modified initial conditions for the 
perfect-fluid evolution [16[ — instead of the initial energy 
density profiles in the transverse plane obtained from the 
Glauber model one should use the Gaussian profiles (with 
the same widths as in the Glauber calculation). 

In this paper we use the standard Glauber initial con- 
ditions and check how our modified dynamics (i.e., the 
dynamics including the transverse-hydrodynamics stage) 
affects the HBT radii. We find that our model offers a 
good description of the HBT radii which are reproduced 
at the level of about 10%. The pion spectra and V2 are 
described with the accuracy of about 10% and 20%, re- 
spectively. This means that our approach provides a very 
good description of the pion production at RHIC. 

The global fit to the RHIC data (including kaons and 
protons) reproduces the soft hadronic observables with 
the accuracy of about 20%. The worst model result are 
obtained for the proton elliptic flow, which is 50% larger 
than the measured value. In our opinion this is caused by 
the neglecting of the hadronic rescattering phase (or the 
effects of the shear and bulk viscosities as shown in p7|). 
In any case, the overall good description of the data ob- 
tained in our model is encouraging, especially from the 
point of view that the problem of early thermalization 
is circumvented and the inconsistencies between the the- 
oretical and experimental HBT results are strongly re- 
duced. 

Throughout the paper we use the three-dimensional 
(3D) and two-dimensional (2D) densities of various phys- 
ical quantities. Thus, to avoid possible confusion we have 
introduced the subscripts "3" or "2". For example, the 
3D energy density (energy per unit volume) is denoted 
by £3, while the 2D energy density (energy per unit area 
in the transverse plane) is denoted by £2. 

The rapidity and spacetime rapidity are defined by the 
expressions 

1 E+ P[l l t + z 

2/=o hl p ' 7 l=o ln 7 ' 1 

2 hi — pit 2 t — z 

which follow from the standard parameterization of the 
four-momentum and spacetime coordinate of a particle, 

p' 1 = (E,p ± ,p\\) = (mxcoshy,p ± ,m_i_sinhy), 

= (t, xx, z) = (t cosh 77, xx, t sinhr;) . (2) 

In Eq. @ the quantity m_i_ is the transverse mass 



m_L = Jm 2 +p 2 x +Py, (3) 

and r is the proper time 

r = v^ 2 - z 2 . (4) 
Everywhere we use the natural units where c = 1 and 

h = 1. 



II. INITIAL CONDITIONS 
A. Initial transverse profiles 

For the boost-invariant systems, which are studied be- 
low, one usually assumes that either the initial entropy 
density, 031 (x_i_) = 03(71, xj_), or the initial energy den- 
sity, £3i(xx) = £3(71, xj_), are directly related to the den- 
sity of sources of particle production, p sr (xi). These 
sources are identified with wounded nucleons or binary 
collisions [l8j . The symmetry with respect to the Lorentz 
boosts along the collision axis means that it is sufficient 
to consider all these quantities in the plane z = 0. 

More generally, a mixed model is used [l9| , with a lin- 
ear combination of the wounded- nucleon density w(xj_) 
and the density of binary collisions n (xi ) . This leads to 
the two popular choices: 

1 — K 

CT3i(xj_) OC /3 sr (xj_) = — - — «j(xi) +Kn(xi) (5) 

or 

1 — K 

£3i(xjJ oc /O sr (xx) = - w(xx) + nn(xx) ■ (6) 

The distributions W(xj_) and 7l(xjJ are calculated for a 
given centrality class from the Glauber model. Follow- 
ing the PHOBOS studies of the centrality dependence of 
the hadron production [20T ] one uses frequently the value 
K = 0.14[2l|. 

By the initial entropy or energy density we understand 
the values of these quantities at a certain early proper 
time r = t, > (the time r = corresponds to the 
moment where the two Lorentz contracted nuclei pass 
through each other). In the standard hydrodynamic ap- 
proaches cr3i(x_|_) and £3i(xj_) characterize the local equi- 
librium properties of the system and r = rj is treated as 
the thermalization time, i.e., the time needed to establish 
local thermodynamic equilibrium. The time evolution for 
r > T\ is fully determined by the hydrodynamic equations 
supplemented with the appropriate equation of state. 

Contrary to standard approaches, in this paper we as- 
sume that the system is not completely thermalized at 
t = 71 . We assume that only transverse degrees of free- 
dom are thermalized at that time, while the longitudinal 
motion is described by free streaming. In this case, we 
may use the following relations connecting the 3D and 
2D densities 

03i(XjJ = — C2i(xj_), £3i(xj_) = — £2i(xj_), (7) 
7i Ti 

where hq is the normalization constant whose physical 
interpretation will be given below. Equations ([5]), (|6l). 
and (O suggest the following form of the initial conditions 
which may be used in our case: 

1 — K 

C72i(xjJ OC — ffl(xi)+d(xi), (8) 
1 — K 

£ 2 i(xj_) OC W(xi) + K)l(xi). (9) 
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The initial transverse flow at r — t\ is set equal to zero. 
The dynamics of the system in the proper time interval 
t; < r < T tr is determined by the equations of transverse 
hydrodynamics. Only at r = r tr a transition to the fully 
equilibrated phase takes place. 



our framework of the transverse hydrodynamics differs 
from its first formulation introduced by Heinz and Wong 
I20L [3011 . These differences are discussed in greater detail 



B. Apparent transverse thermalization 

Before we discuss the concept of transverse hydrody- 
namics, it is useful to recall the arguments for purely 
transverse thermalization of the systems produced in 
hadronic and nuclear collisions. 

In the string models, the production of quarks and glu- 
ons is often understood as the result of the Schwinger 
tunneling mechanism which leads to the Gaussian 
transverse-momentum distributions. However, as shown 
by Bialas, the fluctuations of the string tension change 
the Gaussian dependence into the exponential one 111. 
The exponential transverse-momentum distributions re- 
semble the Boltzmann distributions and the produced 
systems may be interpreted as "transversally" thermal- 
ized. Of course, in this case the thermal character is 
not a consequence of multiple scattering of the produced 
partons but the result of the specific mechanism of parti- 
cle creation (in this context we may talk about apparent 
thermalization) . 

Analogous features appear also naturally in the color- 
flux-tube models, where the color fields oscillate (22l.[23|. 
Moreover, in this case the tunneling particles have no 
longitudinal momentum (in the pair rest frames) (22j, 
which implies that the initial "longitudinal" pressure is 
zero. 

A similar physical picture is present in the theory of 
Color Glass Condensate 24] . In this case, for very early 
proper time r <C 1/Q S) where Q s is the saturation scale, 
the classical gluon fields lead to the following energy- 
momentum tensor 12 5l 12611 



-«i/Q s 



/e 3 (r) 

e 3 (r) 

e 3 (r) 

V 




(10) 



At later proper times, r ^> 1/Q S , both the analytical 
perturbative approaches [27| and the full numerical sim- 
ulations [23 lead to the form 



->l/Q s 



/6 3 (r) 

e 3 (r)/2 

6 3 (r)/2 

V 00 



(11) 



The form of the energy-momentum tensor (jlip and 
the physics of the string models suggest that during the 
early evolution of matter the longitudinal pressure is sig- 
nificantly lower that the transverse pressure. These ob- 
servations trigerred the development of the concept of 
the transverse hydrodynamics [12, [l3j]. We note that 



III. TRANSVERSE HYDRODYNAMICS 

The equations of the transverse hydrodynamics follow 
from the energy-momentum conservation law, 



W = 0, 



(12) 



with the energy-momentum tensor defined by the formula 



If = — [(6 2 + P 2 ) U"U V - P 2 (g^ + VV U ) ] . (13) 

T 

The normalization constant no may be interpreted as the 
density of transverse clusters in rapidity. The clusters are 
formed by groups of partons having the same rapidity. 
They are 2D objects, whose thermodynamic properties 
are described by the 2D thermodynamic variables: e 2l P2, 
02 and T2 (2D energy density, pressure, entropy density, 
and temperature, respectively). These quantities satisfy 
the standard thermodynamic identities: 

£2 + ^2 = T 2 a 2 , 
de 2 = T 2 da 2l 

dP 2 = o 2 dT 2 . (14) 

The baryon chemical potential is neglected here, since we 
consider the midrapidity region where the baryon density 
is very small. The definition of the energy-momentum 
tensor (IT3l) contains the two four-vectors, 

= (uq cosh 7/, u x , u y , uq sinh?7), 
W = (sinhr/, 0,0, cosh??), (15) 



where 



= (u°,u x ,0) = (u°,tix,u»,0) (16) 



is the hydrodynamic flow in the plane z = 0, while 77 is the 
spacetime rapidity defined by Eq. ([1}. The four- vectors 
and are normalized in the following way 



1. 



v% = 



■1, 



WVp = 0. 



(17) 



The four- vector C/ M combines the motion of the fluid el- 
ement in a cluster with the motion of the cluster, thus 
it corresponds to the flow four-velocity in the standard 
hydrodynamics. The term V^V U in (fTB")) is responsi- 
ble for vanishing of the longitudinal pressure, i.e., in 
the local rest-frame of the fluid element, where we have 
U» = (1, 0, 0, 0) and V = (0, 0, 0, 1), one finds 



(£2 0\ 
P 2 
P 2 

V 00/ 



(18) 
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One can notice that exactly the same structure of the 
energy-momentum tensor appears in the theory of the 
color glass condensate and glasma for t>1/Q s , see 
Eq. (EE) 

We solve the equations of the transverse hydrodynam- 
ics numerically using the equation of state 2 



Pi 



£2 
2 



^ 3 C(3) 
2tt 



<7^3 
J 2i 



(19) 



where v g — 16 reflects the spin and color degeneracy of 
gluons. Our investigations are restricted to the midra- 
pidity region (z rs 77 0), where the partonic system 
may be treated as boost-invariant. 

The structure of the energy-momentum tensor (|13|) 
suggests that it is convenient to introduce the three- 
dimensional densities of the transversally thermalized 
system, 



n 



-£2, 



^3 



n 



-cr 2 , 



(20) 



With the help of this notation, the equations of trans- 
verse hydrodynamics take the form 



(<T*U») = 



0. 



(21) 



Here the upper equation is the analog of the Euler equa- 
tion in classical hydrodynamics, while the lower equa- 
tion describes the conservation of the entropy. Using 
the equation of state (|19[) we may rewrite the entropy- 
conservation law in the form 



^p°T|t^) =0. 



(22) 



The structure of Eqs. (|21j) and (|22|) indicates that they 
are scale invariant — the flow profile does not change if 
the temperature is multiplied by an arbitrary constant. 
We shall use this property later to match properly the 
transverse hydrodynamics with the standard perfect-fluid 
hydrodynamics. 



IV. LANDAU MATCHING CONDITIONS 



equilibration is a gradual process. In our effective model, 
however, we treat this transition as a sudden change from 
the transverse to the standard perfect-fluid hydrodynam- 
ics. We stress that this approach is an approximation 
where the continuous process is replaced by a delayed step- 
like change. 

The sudden equilibration transition is described by the 
Landau matching condition 



1 3 u "i 



(23) 



where T^ v is the standard energy-momentum tensor of 
the perfect-fluid hydrodynamics 3 



P39 



(24) 



Here £3 and P3 are the three-dimensional (3D) energy 
density and pressure of the system immediately after the 
equilibration transition. We assume that they are given 
by the equation of state constructed in Ref. [32[ . 
Equations (|23j) and (JUJ) give 



n 



-£■2 — £3, 



(25) 



which should be supplemented by the requirement of the 
entropy growth, 



tr "0 . 

>3 = <?2 < (J3, 

T 



(26) 



where S3 is the 3D entropy density. Dividing both sides 
of Eqs. (|25|) and (|26l) one obtains 



To > 



3 £3 

2<7 3 ' 



(27) 



Certainly, our treatment of the full equilibration tran- 
sition is very much simplified. More elaborate approaches 
would describe this kind of transformation using kinetic 
theory or dissipative hydrodynamics (33l - [35| . For exam- 
ple, in the dissipative hydrodynamics the viscosity cor- 
rections increase transverse pressure and decrease longi- 
tudinal pressure. In the boost-invariant one-dimensional 
case we find 13311 



The large collection of the RHIC data (transverse- 
momentum spectra, ratios of hadronic abundances, etc.) 
suggests that the system produced finally in Au+Au col- 
lisions at the highest beam energy of 200 GeV per nu- 
cleon pair is very well thermalized. For our modeling this 
means that, at a certain stage, the initial transversally 
thermalized system must undergo the full 3D equilibra- 
tion. In a microscopic approach we expect that the 3D 






V 



Pa 



2 2 

3 T 



P-3 



111 

3 r 








P3 - 



\ 



4 r/_ 

3 T 



(28) 



where 7/ is the shear viscosity. The matching of the tensor 
(|28p with the forms (|10p or (fTl"]) requires, however, very 
large values of 77. This observation questions strict ap- 
plicability of the dissipative hydrodynamics at the very 



2 The form of the equation of state used in this paper is 
valid for bosons. For the classical massless particles on 
the plane P2 = v a T% J (27r), and for the massless fermions 
P 2 = 3y 9 T 2 3 C(3)/(87r) [p. 



3 The use of the same four-vector in the definitions of T^ v and 
T^ v ensures that the energy flux on both sides of the transition 
surface has the same direction. 
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early stages. Thus, we are of the opinion that our ap- 
proach based on the use of the transverse hydrodynamics 
offers an interesting alternative for viscous hydrodynam- 
ics in this case. 

As shown below, the important aspect of our proce- 
dure is that the sudden equilibration transition may take 
place at Tt r of about 1 fm. In our opinion, this result in- 
dicates that the "equivalent" gradual equilibration pro- 
cesses may be extended in time and no assumption about 
the sudden and full thermalization of the system at very 
short times (r; of about 0.2 fm for perfect-fluid hydrody- 
namics) must be made. 



V. PERFECT-FLUID HYDRODYNAMICS AND 
FREEZE-OUT 

We assume that after the transition time rt r the mat- 
ter behaves like a perfect fluid. The spacetime evolution 
is described with the help of the formalism developed 
in Ref. [Hj]. Its characteristic feature is the use of the 
modern equation of state which interpolates between the 
lattice QCD results and the hadron-gas calculations [HJ . 
We note that a similar equation of state has been con- 
structed very recently in Ref. (37j . 

In a certain way, we may treat the Landau matching 
conditions as a method to deliver initial conditions for 
the perfect-fluid hydrodynamics at t = r tr . Since the 
system has spent already some time in the expanding 
"transverse phase" , the transverse flow at r = r tr is not 
zero (the issue of the initial non-zero transverse flow is 
discussed in [38l - l40j V In fact, this flow is relatively large 
because the development of the transverse flow in the 
transverse hydrodynamics is fast (larger sound velocity, 
no losses of energy due to the longitudinal work, etc.). 
Of course, the comparison with the data eventually con- 
straints the amount of the transverse flow that is allowed 
to be generated in the transversally thermalized phase. 

The physical observables are obtained by using the 
Cooper-Frye formula with the hypersurface defined 
by the constant value of the freeze-out temperature, 
X3f = const. Thus, we adopt the single freeze-out sce- 
nario and do not distinguish between the chemical and 
thermal frcczc-outs 411. This is of course an approxima- 
tion that is known to work well for pions and worse for 
protons @, [HI ■ We note that only the fully equilibrated 
region of the spacetime is used to construct the freeze-out 
hypersurfce (in other words, there is no hadronic emission 
from the transversally thermalized stage for r < Tt r ). 

The thermal Monte-Carlo code THERMINATOR ^ 
is used to generate primordial particles that include 
stable hadrons and all known hadronic resonances. 
THERMINATOR simulates the decays of resonances which 
proceed in cascades. All studied observables (spectra, 
the elliptic flow v 2 , and the HBT radii) are calculated 
with the help of the Monte-Carlo method. In particular, 
the femtoscopic observables are obtained with the help 
of the two-particle Monte-Carlo method, where the cor- 



relations between the pions are introduced by the appro- 
priate weights defined by the squares of the two-particle 
wave functions 1431. 



VI. RESTRICTING THE MODEL 
PARAMETERS 

Our model has altogether nine parameters. The main 
six parameters are: no - the overall normalization, T 2l - 
the initial central temperature of the transversally ther- 
malized system, - the freeze-out temperature, tj - 
the initial proper time, r tr - the 2D — > 3D transition 
time when the Landau matching conditions are applied, 
and k - the parameter describing the admixture of the 
binary-collision density in the initial conditions (|8|) and 
©. Similarly to Ref. [4J|, we use the freeze-out temper- 
ature 



T 3f = 140 MeV. 



(29) 



In the thermal generation of the particles on the freeze- 
out hypersurface, the extra three chemical potentials are 
also introduced: /is - the baryon chemical potential, us ~ 
the strangeness chemical potential, and u/ 3 - the isospin 
chemical potential. The values of the chemical potentials 
are much smaller than T 3 f , hence, their effect on the evo- 
lution of matter is neglected. They appear only in the 
thermal distribution functions used to generate particles 
on the freeze-out hypersurface. Following Refs. 45[ we 
use the values 

fi B = 28.5 MeV, fi s = 6.9 MeV, m 3 = -0.9 MeV. (30) 

Our numerical calculations were performed for the ini- 
tial conditions of the form ([5]) and ©. A better agree- 
ment with the data (especially for the elliptic flow) was 
obtained generally for the case ©, hence, in the follow- 
ing we shall restrict our considerations to the situation 
where 



7T 

oc -w (x±) + KfT (x±) . (31) 



£2 (Ti,Xj_) 



The normalization constant required in pip deter- 
mines the 2D initial central temperature of the system, 
T 2i = T 2 { n ,0). 

The selection of the optimal values of T 2 i and n in (f3"Tj) 
will be discussed in more detail below. We want to stress 
now the fact that as long as we ignore the Landau match- 
ing condition for the entropy our results are insensitive 
with respect to the following rescaling of the parameters 
n and T 2 i, 



n ^Xn , T 2i ^A- 1/3 T 2i . 



(32) 
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FIG. 1: (Color online) Contour plots of the quantities A describing the average relative differences between the model results 
and the data, see Eqs. (|36[) and (1371) . The values of A are shown as functions of k and T2i. The two upper plots show the 
agreement for the pion and proton spectra, while the two middle plots show the agreement for the pion+kaon and proton vi. 
The lower plot on the left-hand-side shows the agreement for the HBT radii, and the lower plot on the right-hand-side shows the 
agreement for the all considered observables (spectra, V2, and the HBT radii). The calculations were done with the parameters 
specified by Eqs. (|29p . (|30|) . (|34p , and 1)351) for the centrality class c = 20 — 30% (the impact parameter b — 7.16 fm). The data 
used to calculate A were taken from Refs. I(i i IK. 



This transformation does not change the energy density other hand, it affects the entropy since we have 
and flow of the transversally thermalized system. On the 

^a 2 ^\V*^a 2 . (33) 
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Thus, if a satisfactory description of the data is achieved 
with certain values of n and T 2 j (with the conserved 
energy during the 2D — > 3D transition) we may always 
use the transformation (|32l) to change the entropy in the 
purely transverse stage in such a way that it grows dur- 
ing the equilibration transition. By decreasing A we ef- 
fectively reduce the number of clusters and make them 
hotter. The decays of such "more organized" and hotter 
clusters produce more entropy. 

Following our previous studies 0, [TH where the trans- 
verse hydrodynamics was followed by the sudden equili- 
bration and freeze-out (without the perfect-fluid stage) 
we used the value 

n = 0.43, (34) 

which turned out to be satisfactory also in the present 
studies. 

We have also found that the optimal values of the pa- 
rameters Tj and Tt r are 

Ti = 0.25fm, T tr = 1.0fm. (35) 

The evolution times of the transversally thermalized sys- 
tem that are longer than 0.75 fm lead to too strong radial 
flow. On the other hand, a slightly better description 
of the data is achieved if the transverse-hydrodynamics 
stage is implemented earlier. Such constraints lead to the 
choice (l35|) . We note that the time interval of the purely 
transverse hydrodynamics considered in this paper is the 
same as the time for the initial free streaming considered 
in Ref. 0]. 

The quality of the fits discussed in this Section is char- 
acterized by the parameters A defined as the average rel- 
ative difference between the model results and the data 

1 N l-th cxpi 
i—l 1 

Here a^ xp is the ith experimental point (e.g., the value 
of the transverse-momentum spectrum or the value of 
the elliptic flow) and xf 1 is the corresponding theoretical 
value. The quantity N is the number of the experimental 
points (e.g., the number of the points in the measured 
transverse-momentum spectrum). The values of A are 
calculated for different observables separately, and also 
for combined observables. In the latter case we define A 
as the average 

A bsl+obs2 = ^ (Aobsl + A Q bs2) • (37) 

We stress that the lack of knowledge about the system- 
atic errors (in the case of the transverse-momentum data 
presented in (46|) or discrapancies between the system- 
atic errors given by different experiments (e.g., the HBT 
results given in (48| and |49| makes a thorough \ 2 analysis 
quite difficult. It requires making comparisons between 
different experiments and such a study goes beyond the 
problems discussed in the present paper. 




r [fm] 



FIG. 2: (Color online) The 2D and 3D isotherms for <f> = 
(upper part) and <f) = it/2 (lower part). The values of the 
parameters are specified by Eqs. ([29}, (J30j , ([34} . ([35l) . and 
d}. The centrality class c = 20 - 30%. 
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r [fm] 



FIG. 3: (Color online) The energy density at the transition 
time Ttr plotted as a function of r for <f> — (dashed line) and 
4> = 90° (solid line). The parameters the same as in Fig. [2] 



A. Determination of T2i and k 

In order to illustrate the choice of the optimal values of 
Ti i and k, in Fig.[T]we show the values of A calculated for 
different physical observables. As discussed above, the 
values of other parameters were motivated by our earlier 
hydrodynamic studies and are given by Eqs. ([29}, (|30[) . 
(|34[). and ([33]) 4 The model calculations are compared 



4 The multidimensional fits using all parameters are very much 
time consuming and not possible for us. In particular, the cal- 
culations of the HBT radii are very long. 
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FIG. 4: (Color online) The entropy density just before (solid 
lines) and after (dashed lines) the transition time r tr plotted 
as a function of r for <f> = and <j> = 90°. The parameters the 
same as in Fig. [2] 



with the data published in Rcfs. [46-48]. We consider 
the centrality class c = 20 — 30% (the impact parameter 
b = 7.16 fm) 

The upper two plots of Fig. [T] show the quality of fit- 
ting the transverse momentum spectra. The experimen- 
tal PHENIX results used to calculate A arc taken from 
Ref. [46l |. We find that the pion spectra are described 
with the accuracy of about 10-60% depending on the val- 
ues of Ti \ and k. A similar accuracy may be achieved in 
the fits of the proton spectra, where A = 10-40%. The 
interesting point is, however, that the pion and proton 
spectra favor different values of T21 and n. It turns out 
that it is impossible to achieve the accuracy better than 
20% for both pions and protons. 

The middle two plots show the quality of fitting the 
elliptic flow coefficient V2- The PHENIX data are taken 
now from Ref. [471] . In the considered ranges of the pa- 
rameters T21 and k we find that the pion+kaon V2 is 
reproduced at the level of about 10-30%, while the pro- 
ton V2 is reproduced much worse, at the level of about 
30-60%. We note that such high values of the proton 
elliptic flow are characteristic for the calculations using 
the crossover phase transition [2l], |37j]- The difficulties 
connected with the correct reproduction of the proton V2 
are illustrated in the panels (a) and (d) of Fig. Q] where 
one can see that the good fits of the pion spectra and 
proton V2 are anticorrelated. 

In the lower left part of Fig. Q] we show the agreement 
between the model and experimental results for the HBT 
radii. In this case the STAR data from Ref. [48[ are used. 
Generally, we find a very good description, at the level of 
10% for all the values of T 2 i and k. The good description 
of the HBT results is a consequence of many factors [5(| ■ 
In our case they include the use of the semi-hard realis- 
tic equation of state [32[ and the implementation of the 
"transverse" phase that speeds up the time evolution of 
the system. 

The analysis presented in this Section suggests that 




012345678 
r [fm] 



FIG. 5: (Color online) The transverse flow at the transition 
time Ttr plotted as a function of r for <j> — (dashed line) and 
<p = 90° (solid line). The parameters the same as in Fig. [2] 



the optimal values of the parameters T 2 i and k are 

T 2i = 530 MeV, 
k = 0.25. (38) 

These two values together with the earlier specified pa- 
rameters, see Eqs. ([2"9"|). (|3TJ)) . and will be used 
in the following presentation. 

We note that Fig. [T] illustrates quantitatively difficul- 
ties of obtaining a consistent good fit to all studied ob- 
servalbles — different physical quantities are described 
well by different ranges of the parameters. 



VII. SPACETIME BEHAVIOR OF 2D AND 3D 
THERMODYNAMIC QUANTITIES 

Figure [5] shows an example of the spacetime evolution 
of matter in the vicinity of the transition time r tr . The 
values of the parameters are specified by Eqs. (j2"9"1) . (jBTJI) . 

([331) . and ([3"8"f. We recall that they correspond to 
the centrality class c = 20 — 30% (the impact parameter 
b = 7.16 fm). In this case the transverse hydrodynamics 
starts at t\ = 0.25 fm and continues till r = r tr = 1 fm, 
when the transition to 3D evolution takes place. 

The upper and lower parts of Fig. [2] show the 2D and 
3D isotherms for = and 4> = 90° > respectively. The 
energy density which is the same on both sides of the 
transition is plotted in Fig. [3] 

The application of the Landau matching conditions im- 
plies that there is no simple connection between the 2D 
and 3D isotherms at r — r tr . However, if the hadronic 
gas was described by the equation of state of mass- 
less particles then Eq. (|27|) would yield the condition 
T 2 > (9/8)T 3 . Thus, we expect qualitatively that T 2 
should be larger than T3 in the transition region. By in- 
spection of Fig. [2] we see that this condition is fulfilled 
for the central part of the fireball (for r < 5 fm if = 
and for r < 6 fm if <fi — 90°). 
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FIG. 6: (Color online) The comparison of the model and ex- 
perimental transverse- momentum spectra of pions, kaons, and 
protons. The data are taken from Ref. [4g|. The model pa- 
rameters are defined by Eqs. (gSJ, fl3DJ, (EH), and (|gg|). 
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FIG. 7: (Color online) The comparison of the model and ex- 
perimental elliptic flow coefficient vi for pions+kaons and pro- 
tons. The data are taken from Ref. [43]. The model param- 
eters the same as in Fig. [6] The bands indicate the error of 
the Monte-Carlo calculations. 



The exact jump of the entropy density is illustrated in 
Fig. SI where the solid lines describe the entropy den- 
sity just before the 2D — > 3D transition and the dashed 
lines describe the entropy just after the transition. In 
agreement with the expectations formulated above, we 
observe that the entropy increases substantially in the 
central part of the fireball. In the outer parts of the fire- 
ball the entropy densities in the transverse and perfect- 
fluid phases are very much similar. For large values of r 
the entropy in the transverse phase becomes smaller and 
the Landau matching condition is violated. However, 
this region is excluded from the physical considerations 
as it corresponds to the 3D temperatures smaller than 
the freeze-out temperature. We note that the entropy 
in the transverse phase can be made even smaller if the 
scaling (|32|) with A < 1 is performed. 

For completeness, in Fig. [3] we present the flow profiles 
at the transition time. By construction the flow is the 
same on both sides of the transition surface. 
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FIG. 8: (Color online) The comparison of the model and ex- 
perimental results for the pionic HBT radii. The data are 
taken from Ref. [48! ]. The model parameters the same as in 
Fig. [6] 



VIII. COMPARISON WITH DATA 

The methods presented in the previous Sections de- 
livered us the optimal parameters for fitting the RHIC 
Au+Au data for the centrality class 20-30%. In this 
Section we present the direct comparison of our model 
results with the data. In Fig. [6] we show the model 
and experimental transverse-momentum spectra of pions, 
kaons, and protons. The averag e ag reement between the 
model and the PHENIX data [46| is quite satisfactory 
(A = 13% for pions, 7% for kaons, and 28% for protons). 

Figure [7] shows our results for the elliptic flow com- 
pared with the PHENIX experimental results from Ref. 
[47| . We are able to reproduce very well the experimental 
pion+kaon elliptic flow (at the level of 12%), however the 
model proton w 2 is too large by about 50%. The large 
theoretical values of the proton elliptic flow are connected 
with the lack of the hadronic interactions in the final state 
or with the neglecting the viscous effects (as pointed out 
recently in Ref. [17j). 

In Fig. |S] we show our results describing the HBT 
radii and compare them with the STAR results [48J . We 
observe a very good agreement between the model and 
the experimental data. The radii are reproduced with 
the accuracy better than 10%. 

We note that a similar analysis may be performed for 
the central collisions. With a slightly changed values of 
the input parameters we are able to reproduce the ex- 
perimental data with a similar accuracy as in the case of 
the non-central collisions. 
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IX. DISCUSSION AND CONCLUSIONS 

The main conclusion of our analysis is that it is possible 
to obtain a satisfactory description of the soft hadronic 
RHIC data in a hydrodynamic model that avoids difficul- 
ties connected with the concept of very fast 3D thermal- 
ization. In our model, the initial evolution of matter is 
described with the help of the transverse hydrodynamics 
that uses the assumption of the transverse thermaliza- 
tion only. The initial thermalization of the transverse de- 
grees of freedom may be naturally explained in the string 
or color-flux-tube models. By combining the transverse 
hydrodynamics with the standard perfect-fluid hydrody- 
namics we include in the schematic way the effects of full 
3D thcrmalization. 

We emphasize that our approach goes beyond the ap- 
plications of viscous hydrodynamics, since the large ini- 
tial anisotropy of pressure may be explained only by in- 
troducing a very large shear viscosity, see Eq. (I28|) . which 
makes the very concept of the kinetic coefficients ques- 
tionable. 

Our model reproduces well the HBT radii indicating 
that the two main problems of the soft physics at RHIC 
may be circumvented. The remaining issue of the large el- 
liptic flow of protons may be very likely explained by the 
inclusion of the shear and bulk viscosity in the standard 
hydrodynamic stage or by the inclusion of the hadronic 
rescattering at freeze-out. The results for other centrali- 
ties will be published elsewhere [5lj . 

In our earlier work [3, EE] we considered a physical sce- 



nario where the transverse hydrodynamics was followed 
by a sudden isotropization and freeze-out. The main dif- 
ference between the approach used in [13, EBl and the 
present model is that we include now an extended phase 
governed by the equations of perfect-fluid hydrodynam- 
ics. The results of the two approaches are, however, 
quite similar (the approach with sudden isotropization 
and freeze-out gives slightly smaller HBT radii, on the 
other hand, their kx dependence is reproduced very well 
yielding a very good description of the ratio -R ut/-Rsidc)- 
The similarities between our present results and the re- 
sults of Refs. [13, EH indicate that the details of model- 
ing the 2D — > 3D transition are not very much important. 
The important point is that the matter should thermalize 
before freeze-out. Of course, the more realistic 2D — > 3D 
switching may take place under different circumstances 
than those discussed in this paper or in Refs. fli]. fl5|. 

The present formalism may be generalized in many 
ways. In particular, the perfect-fluid stage may be re- 
placed by the viscous hydrodynamics. The inclusion of 
the viscous effects is also possible in the transverse hy- 
drodynamics. This opens new possibilities for the appli- 
cations of the concept of transverse thermalization and 
for studies of the effects of anisotropic early dynamics. 
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